SOME PARALLEL ALGORITHMS IN COMPUTATIONAL GEOMETRY 


A Thesis Submitted 

m Partial ulfilment of the Requirements 
for the Degree of 
MASTER OF TECHNOLOGY 


By 

V PADMAVATHI 


to the 

DEPARTMENT OF COMPUTER SCIENCE AND ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

FEBRUARY 1991 



CERTIFICATE 



This IS to certify that the Thesis worl- entitled Some 
Parallel algorithms in Computational Georretry has been carried 
out under rry supervision and has not been subrritted elsewhere 
for the award of a degree 


Station KANPUR 
Date 4 th Feb 91 


(Dr Asish Mul- hopadhyay i 
Assistant Professor 
Dept of C S E 
I I T Kanpur 


(llj 



n Q ^ r‘r> f QQ j 

entral library 

I I T KANPUR 

4cc No A IJ Qf\8t) 

C<;e iW' ^ PAP-SotM 




CERTIFICATE 



This IS to certify that the Thesis work entitled Some 
Parallel algorithms in Computational Georretry has been carried 
out under rry supervision and has not been submtted elsewhere 
for the award of a degree 

V "4'“ V ^ 

Station KANPUR < Dr Asish Nuthopadhyay 

Date 4 th Feb 91 Assistant Professor 

Dept of C S E 
I I T Kanpur 


(iU 



abstract 


Computational 6eo-.atry addreasas algorithmic problems in 
diverse areas It is seldom obvious hou to generate parallel 
algorithms in this area since most popular techniques involve an 

essentially sequential approach 

Here parallel algorithms have been proposed for certain 

problems in computational Beometry using CREMPRAH processors 
Among these are the linear programming problem in two dimensions 
the minimum spanning circle problem for a set of points in the 
Plane for a set of line m the Plane and also for a set of 


1 1 4.^ pjav;:illsl algorith(TS also have been 

planar line segirents paraiiei di-i 

p XL. 4 4 wi +hp cihortest line segrrent frorr 

for the corrputation of tne snurt.c=>i. 

given conve polygon is weatly eternally visible and 
reconstruction of an orthogonal polygon given a circula 
of its visibility graph 


designed 
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for the 
err bedd a ng 
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CHAPTER 1 


INTRODUCTION 

1 1 CotTputat lonal Geoiretry and Parallel Algor ithtrs 

CoiTputat lonal Geoiretry addresses algorithiric problerrs in 
diverse areas A large nurrber of application areas such as pattern 
recognition corrputer graphics irrage processing operations 
research statistics corrputer aided design robotics etc have 
been the incubation bed of this discipline since they provide 
inherently georretric problerrs for which efficient algorithms have 
to be developed This field has been flourishing because of the 
rich interaction with many application areas 

Visibility IS one of the most fundamental topics in 
Computational Ceorretry Visibility problems find application in 
irany areas such as graphics automated cartography image 
processing and robotics Also visibility problems often appear as 
subproblems of many other problems in computational geometry <like 
the shortest path problems with obstacles^ 

Since 1975 there has been a rapid development of sequential 
algorithms for geometric problems but until 1985 there was little 
published about parallel algorithms for such problems 

It IS seldom obvious how to generate parallel algorithms in 
this area since popular techniques include an eiplicitly 
sequential (iterative' approach 



Here parallel algorithrrs have been developed for probleiTB in 
visibility and opt ifr isat ion 

1 2 Survey 
Serial Algorithms 

There has been a wide developrrent of sequential algorithrrs 
for georretric problerrs We report here a few results that are 
particularly relevant to our worf 

The notion of weah visibility has received attention in both 
tr atherr at 1 cs and corrputer science literature Avis and Toussaint 
CAT813 showed that given a sirrple polygon P and a specified edge 
e whether P is edge visible frorr e can be deterrr ined in Dtn^ 
tirre Sact and Sun CSS891 discovered a linear tiire algorithrr for 
deterrTining all such edges of a given polygon Bhattacharya CBTQ93 
gave a linear tirre algorithrr for corrputing the shortest line 
segrrent frorr which a given conve polygon is wea^ly externally 
visible 

A linear tirre algorithir has been given in CQJ87I1 for the 
reconstruction of an orthogonal polygon given a circular errbedding 
of its visibility graph 

Nirrrod Meggido CMN83] gave a linear tirre algorithrr for linear 
2 3 

prograrrir ing in R and R 
Parallel Algorithrrs 

Not rruch has been published on parallel algorithrrs for 
georretric problens Sorre of the results in parallel algorithrrs in 
Corrputat lonal Georretry are discussed below 



Atallah Cole and Goodrich showed that one can corrpute ttie 
portion of the plane visible frorr a point of the plane in the 
presence of a collection of n line segrrents in OClog n^ titre using 
0(rO processors The nurrber of processors needed for this problerr 
was irrproved by Atallah and Chen CACS93 to O(n/log n) for the case 
when the line segrrents forrr a sirrple polygon ElGindy and Goodrich 
CEG873 showed that one can determne the shortest path between two 
points in a sirrple polygon in OClog n 1 tirre using OCn) processors 
and one can construct the shortest path tree frorr a verte v (the 
union of all shortest paths frorr v to other vertices) in OClog n) 
tirre using OCn) processors Aggarwal Chagelle Guibas et 
al CACSOl presented efficient parallel algorithrrs for basic 
problerrs in Corrputat lonal Georretry vi convex hulls voronoi 
diagrarrs detecting line segrrent intersections triangulating 
sirrple polygons minimsing a c i rcurrscr ib i ng triangle in OClog n) 
OClog n) OClog n) OClog n) OClog n) tirre with OCn) OCn) OCn 

log n) OCn log n) OCn) processors respectively Recently 
Goodrich ShaucI- Surranta Guha C6S9(fl showed that a data structure 
called the stratified decorrposit ion tree can be built in OClog n) 
tirre using OCn> processors that allows one to construct an 
irrplicit repr esentat ion of the shortest path between two points 
inside P in OClog n) tirre using OCn) processors It is also shown 
that the visibility graph flfoiii, the vertices of P can be corrputed 
in OClog n) time using OCn log n + ^/log n) processors where [ is 
the number of edges in the graph 

For all these algorithrrs the trodel of cotrputation is the 



CREW PRAM iTodel (Concurrent Read Exclusive Write Parallel RAM ^ 


1 3 Model of Corrputation 

An BIMD (Single Instruction Multiple Data> CREW (Concurrent 
Read E elusive Write) PRAM (Parallel RAM^ rrodel of corrputation is 
used An SIMD corrputer consists of a nurrber of processors 
operating under the control of a single instruction strearr issued 
by a central control unit The processors each have a snail 
private rrerrory for storing prograrrs and data and operating 
synchronously during a given tirre unit a selected nunber of 
processors are active and e ecute the saire instruction each on a 
different data set the rerraining processors are inactive The 
processors are synchronised in the sense that if a set of 
instructions is e ecuted in parallel then each trust be allowed to 
finish before the ne>t set of instructions is started 

In this rrodel there is an unbounded global rrerrory which is 
shared by all the processors Processors comrr un irat e with each 
other through this cotrrron rrerrory No two processors are allowed to 
write on the sarre location of the global rrerrory si rr ul taneousl y 
But any nurrber of proces_.Drs can sirrultaneously read frorr the sarre 
location of the global rrerrory 


Execution of Parallel Algorithtrs A parallel algorithrr is one 
that IS designed to run on a parallel corrputer 


Paral lei 

r unn mg 

t irre 

IS 

the nost irrportant i 

treasure 

in 

evaluating a 

parallel 

algor i thrr 

It IS defined as 

the 

t lire 

required to solve a problerr 

that 

IS the tirre elapsed 

frorr 

the 



rroirent the algorithrr starts to the rrorrent the algorithfr 
teriTinates For a problerr of si e n the parallel worst case 
running tirre of an algorithrr a function of n is denoted by t(nl 

A good indication of the quality of a parallel algorithir for 
sorre problerr is the speedup it produces 

Worst case running tirre of the fastest 1 nown 
sequential algorithrr for the problerr 

speedup = — — 

Worst case running tirre of the parallel algorithrr 

The nurrber of processors required to solve a problen is 
another criterion for assessing the value of a parallel algorithrr 
The nuirber of processors required by an algorithm a function of 
n IS denoted by ptnl 

The cost of a parallel algorithrr is defined as the product of 
the parallel running tirre and the nurrber of processors used For a 
problerr of size n the cost of the parallel algorithn a function 
of n will be denoted by cCnl 
c(n> = p(n' t(n^ 

The efficiency of a parallel algorithrr is defined as 

Worst case running tiire of fastest 1 n awn 
sequential algorithrr for the problerr 

efficiency = — — — — 

Cost of the parallel algorithrr 

One of the rrajor tashs of parallel algorithir design is to 
cone up with parallel algorithrrs that are optirral i e that run as 
fast as theoretically possible for the problerr they solve and 
sirr ultaneously have a tirre ^ processors bound that is within a 
constant factor of the tirre coirple ity of the best Inown 



sequential algorithrr foi' the problerr 


1 4 Overview of the thesis 

In chapter 2 parallel algorithrrs have been proposed for the 
Linear Prograrrrr ing Problerr and related prablerrs 

In chapter 3 a parallel algorithrr is given for the 
corrputation of the shortest line segirent frorr which a given conve 
polygon is wea^ ly e ternally visible 

In chapter 4 a parallel algorithn has been presented for the 
recon truction of an orthogonal polygon given the circular 
errbedding of its visibility trees 

We conclude in chapter 5 surrrrar ising the results obtained 
and indicating directions in which further research can be carried 
out 
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CHAPTER 2 


LPP AND RELATED PROBLEMS 


2 1 Introduction 

Several ^ nown problerrs in corrputational georretry are closely 
related to the problem of linear progratrrr ing The probleir of 
finding the srrallest circle enclosing n given points in the plane 
the problerr of corrputing the rrinirrum spanning circle of a set of n 
given lines line segrrents are all related to the problerr of 
linear progr arrrr ing 

Nirrrod Megiddo proposed a linear tirre algorithrr for linear 
2 3 

pruyrarrriing in and R CMNBd] He also proposed linear tirre 
algorithirs for the rrinirrurr spanning circle problerr of a set of 
points Very recently linear tirre algorithns for finding the 
mnirruiT spanning circle of a set of lines and line segrrents have 
been proposed by Mulhopadhyay Bhattacharya et al CBM9£.i] and 
CJM9ff3 

In section 2 2 we give a parallel algorithir for the Linear 
Prograrrming Problem in two dirrensions In section 2 3 the rrinirrum 
spanning circle problerr for a set of points is solved In sections 
2 4 and 2 5 the problerrs of corrputing the mniiruir spanning circle 
for a set of lines and line segrrents are dealt All the algorithirs 
run in OCn^) tirre on 0(n^ CREW PRAM processors 
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2 2 The Linear Pr ogr arrfr ing Problerr 
2 2 1 The Problerr 

The linear progr arr rr mg problerr in two dirrensions can be 

stated as 

mn c^ ^ + C 2>2 


such 

that a . > . + a _x„ 
111 IZ z 

^ 0^ (1 

= 1 

n ) 


It IS 

convenient for us 

to deal with 

the problerr 

in an 

equivalent 

f orrr 





rr in 

y 





(> y ) 






such 

that y > a + b 

1 1 

(1 € I^ 1 





y < a + b 

X 1 

<1 € I^^ 

a 

< < b 


where 1 1 1 

+ |I^| < n 00 < 

a b < 00 




Define the 

f unct ions 





g( ) 

= rra ^ ^i 

1 € I^} 




h( ) 

= rrin £ a^x + b^ 

1 e I 




Both 

these functions are 

piece wise 

linear and g is 

can ve 

while h IS 

concave A nurrber 

X a S 5 b 

xs 

said to be 

f e asib Is 


if g( ) < h(>^ The problerr can be posed in a one dirrensional 
f or rr 

rrinimse g( ) 
such that gf ) ^ 

a < X S b 

Testing a value of > Given any value x of > (a x < b^ we 
test (i^ Is y feasible (ii^ If x is not feasible then if there 
are any feasible values of s then they rrust all lie on one side of 
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X 


our test either deterrrines that side or concludes that no 


feasible value exists ^iii^ If x is feasible then our test will 
recognise whether is also optirral and if not at will tell us on 


what side 

of > 

the 

rr xnirr urr 

lies 

( see 

fig 2 

2 

Define 

s 

g 

= FT xn 

la 

1 

1 

4E 


a y 

1 

+ 

b 

1 

= g ( X 

^ } 

s 

g 

= ftax 

{a 

1 

1 

€ 

^1 

a > 

1 

4- 

b 

1 

= g (> 

^ } 

^h 

= rr 1 n 

£a 

X 

1 

€ 

^2 

a 

X 

+ 

b 

1 

= h ( X 

> 5 

^h 

= nay 


1 

€ 

^2 

a > 

1 

4- 

b 

1 

= h(x 

^ } 


Case ^ g ( > ) h C 1 

Consider the function f(>^ = g( 1 hf > 

This function is conve so all the values of x such that f< ) 

< 0 lie on one side of In order to tell the correct side we 

lool at the one sided derivative of f at 

If s S, then fC ^ is ascending at > so that a feasible 
g h 

can only be srraller than If S s. then ff 1 is descending at 

g n 

> so that a feasible > can only be larger than The rerraining 

ca^e IS wher s S. < 0 S B s. In this case f attains its 

g h g h 

rrininurr at > i e there are no feasible values of 
Case ii g( ^<'h(>) 

If s '^0 then the optirral solution (denoted by rr ust 

g 

satisfy X* X If S < 0 then x^ x Otherwise s < 0 < S and 

g g g 

X itself IS a mnirrun of g 

Case III gti ) = h(> > 

If s 0 and s ^ S. then \ > 
g g h 
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5a 


If S 0 and S ^ s. then > 

9 9 h 

Otl-erwise > itself is a rrinirrun pf gC>^ under the constraint 
g( ) S h(> ) 

2 2 2 The Algorithm 

s 

The procedure par-lpp solves the LPP problerr in OCn ^ tirre 
1 s 

using 0(n ) CREW PRAM processors Pctr-s&l is used by pax—lpp 

th 

The procedure pcuc—s&l gives the ^ elerrent of the input sequence 

E 1 G 

of length n in OCn ) tirre using QCn ) processors CA'^8‘^11 
pax—s&lCS M-> 

(1) If |B| 3 then using one processor return the f elenent 

e Ise 

1. G 

subdivide S into |S1 subsequences each of 

G 

SI p> |S1 and assign a subsequence to each 

processor 

1 “E 

(2^ For 1 = 1 to |S1 do in parallel 

P^ finds the rredian rr ^ of its associated subsequence using the 

"t h 

sequential procedure select writes rr ^ in MCil the i 

position of the array M in shared rrerrory 

(3) Find the rredian rr of M by calling par-selCM 

(4) Subdivide S into three subsequences S. S^ of eleirents 

X jL. O 

srraller than equal to and larger than rr respectively 

(5) If |Sj| > h then par-selCS^ /O 

if (S. I +-(S^| S: P then return rr 

1 iL 

else par—selCS^ k. ~ jS^ | - jS^p 


else 


end 1 f 


end 1 f 
Analysis 

Step 1 ta^es constant tirre 

6 

step 2 needs Cj^ri tirre (where c is a constant > as the 
sequertial select is of linear tirre 

Step 3 tal-es t(n^ tirre as 1M| = ^ 

Subdividing S into S. S„ S„ can be done by letting each 

12 3 

processor split its associated ssquapca into three lists s/ 
S,,^ and of elerrerts srraller than equal to and larger than rr 

respectively in Q(n^) tirre Then step 4 tales 0(n^l tirre The 

and S«^ are rrerged to forn Si Br, S^ respectively as follows 

iC. C> 1 O 

To rrerge denote For each i ^ = XJs^ (j = 1 to n) 

1 ■*“£ 

IS found by parallel prefi corrputation in O(log n ) tirre 

Tal mg = L) all the processors s irr ul taneousl y write their lists 

in Sj with starting to copy ^ in position i i ^ 

E 

array Hence the tirre required by step 4 is domnated by c^n 

where c,-, is a constant 

i “**G 

Since IT IS the rredian of M n /2 elerrents of S are 
guaranteed to be larger than it Every elerrent of M is srraller 
than n^/2 elerrents of S Thus £ 3n/4 Sirrilarly jS^I^ 3n/4 

Hence step 5 tales t(3n/41 tirre 

t(nl - c^n^ + t(n^ + c^n^ + t(3n/4l = 0(n®l 

The following procedure solves the LPP problerr The problerr 
IS talen to be in the forrr 
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rr 1 n y 

( y^ 

such that ySra^+b fie I,) 

11 1 

y<ax+b fie a < > < b 

1 1 2 

where 1I.|+|I^| in oo<a b<CD 

Arrays A B contain s and s respectively 
poii — IppCI ^ I ^ A B a bS> 

C 1 ^ I "T 1 1 ^ I and 112 I are con tant use the sequential algorithm 

e 

else divide the set I into subsets of si e n each and assign 
a processor to each Simlarly divide the set into subsets 

g 

of 51 e n and as ign a processor to each 

For 1 = 1 to n^ ® do in parallel 

arranges the elerrents in its associated sequence in 

di joint pairs 

X = { } 

1 

For each pair Cl in its associated sequence does 

If a^^ = a^then if <1 j'> e then drop one of the 

constraints 


y > a , + b and 

J J 

y 2: a^ + b^ 

else drop one of the 
constraints 



1 



el«='e add y 


, to the set X 
^ J 1 

<3i Let X = union of all X s 

1 

= par s&lCX \X\/'a:> 

Let C = { } H = { } 

1 6 

For X = 1 to n do in parallel 

P firds f , = a, y + b, for each i- c the associated 
subsequence of 

If the associated subsequence of P^ e then P^ adds to 

the et G 
else adds to H 

(5^ g = pars&lCG \G\^ 
h = par-selCH fJ> 

(6) Let SI = { } 

S = { } 

1 G 

For 1 = 1 to n do in parallel 

For each elerrent ^ in its associated subsequence does the 

following 


If 

a, > + b. 

1 IT t 

~ g then 

add a^ 

to 

B1 

if 

a, > + b. 

V rr [ 

= h then 

add aj^ 

to 

S2 


(71 s 

g 

” pczr-^^lCSl 

i> 

^h 

- pctr-s&lCS2 


s 

g 

= pcLTS&lCSi 

|Sf p 

^h 

- pcxrs&lCSB 

|S^p 


(81 Pj^ checl-5 g and h as fallows 
case 


13 



g 


h 


g h 

if s 

g 

> pi then Ca 

b3 

= Ca 3 

rr 




If S 

g 

< 0 then Ca 

bl 

= C> b3 

rr 




else 

return > 

rr 





g = h 


^ 0 and s 5: 

g 

^h 

then Ca b3 

= Ca 



if S 

g 

0 and S < 
g 

^h 

then Ca b3 

= C> 

n 

b3 


else return > 

rr 


g 

h 1 f s 

g 

S, then 
h 

Ca 

b3 = Ca 

3 

rr 


If S 

g 

s. then 
h 

Ca 

b3 = C 

rr 

b3 


else return infeasible 

1 B 

(9^ For 1 = 1 to n do in parallel 

For each elerrent > e X P tests if > e Ca b3 

1 1 

if « [a bl then drop oite of the constraints 
end for 
end for 

Let I, be the sets obtained after dropping the 

const r a 1 nt 5 

(irn par-lppCI^ A B a b:> 

223 Analysis 

Steps 2 4 S and 9 require 0(n^^ tiire as each processor 

spends a linear arrount of tirre on its input 

Steps 3 5 7 tale O(n^) tifre since this is the tirre talen by 

par-s&l 

Here the forrration of the sets X H Q I 1-^ can be done 
in the lines simlar to the one used to font in par-s&l (as 

described in the analysis of step 4 of pax— sol) and henre can be 
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dcr e in O(n^) tirre 

At least half of the critical values will not be in the 

1 J 

interior of the new interval We will thus be able to drop one 
constraint per each such pair which is at least a quarter of the 
set of constraints f including those that have been dropped prior 
to the evaluation of x ) We are thus left with a LPP in the plane 

IT 

with atrrost 3n/4 constraints Hence step 10 tales tt3n/4l tine 
since after eliirination 3n/4 constraints remain 
Rest of the steps tale constant tine 
t(n) = t(3n/4l + c^n® = 0(n®l 

2 3 Minirrum spanning circle for a set of points 
2 3 1 The ProbleiT 

Given a set of n points (a b 1 i == 1 to n in the Euclidean 

11 

plane the problerr is to find the snallest circle enclo_>ing these 
points Forrrally we are lool ing for a point ( y) so as to 

IT inirr ise 

rra i (tx a )^ + (y b 1 < i < nJ 

1 1 

Thus the point t yi is an optinal location of a facility 
if we wish to iTinirrise the largest distance that a custorrer would 
have to travel fron his residence to the facility 

We first develop a parallel algorithrr for the constrained 
case where the center is constrained to lie on the > a is 
Consider the probleir of trininising 

g(>l = na> { (x - ^ 1 < i S n } 
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2 3 2 The Algorithir 

Lei denote the irinitriscr of g<.K> 

The procedure par -cons -point determnes the center and the radius 
of the rrinirrurr spanning circle (constrained case) of the given n 
points using 0(n^ CREW PRAM processors in O(n^) tire 
par— cons— point CS^ 

(1) If |S1 IS not a constant subdivide B into subsequences of si e 
€ 

n and assign a processor to each subsequence 


(2) X = t } 


For 1 = 1 to n 


1 e 


do in parallel 


pairs up the points in its associated subsequence 


For each pair (a 


a in its associated subsequence 

J + 1 


if a 


a , , drop one of the points (a b > and 
j + 1 J J 


(a 


J+1 


j+ 1 


coiTPute > J 
add > , ^ to X 

end for 
end for 


(a 


J + 1 


a +b 
J J + 1 


b / 
J 


: ( a 


j + 1 


(3) = par-s&lCX 

F = { } 


P does 
1 


a ) 
J 


£ 

For 1 = 1 to n do in parallel 

2 2 

P corrputes f = (x a ) + b for each j in its 

1 J IT J J 

associated subsequence and adds it to F 
end for 

(5) g nr = par-selCF |F|^ 
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(6) I 


= { } 

1 G 

For 1 = 1 to n do in parallel 

For each point b.l in its a sociated subsequence P does 

J J ^ 

2 2 

if (> a ) + b = g IT add j to I 

w w 

end for 
end for 

(7) Initialise right = false left = false center = false 

G 

If jl I n use one processor 


else 

€ 

Subdivide I into subsequences of si e n and assign ca h 
of these subsequences to a processor 
For 1 = 1 to n^ ® do in parallel 


P initialises right = false left 
1 1 1 

false 


false centre 


sets right = true if > a for every j e its associated 

1 IT J 

subsequence 


sets left = true if x > a for every j c its a'-socialed 

1 IT J 

subsequence 


end for 

(81 Bet left = true if all left^ = true 

Set right = true if all right = true 

1 

else return (> 0) and g m 

rr 

1 G 

(91 For 1 - 1 to n do in parallel 

P does 
1 

if left = true then if discard (a b > 

J j + l tr J J 
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if right 


= true then if x, ... y discard <a b .. > 

J j + 1 rr j + i j + l 

end for 

S - set of rerraining points 
par-cons-pointCS J> 

Analysis 

Step 1 ta^es constant tiire 

Steps 2 4 6 7 and 9 ta^ e 0(n^' tirre as in each of thf e 

each processor spends a linear arrount of tiire on its input 
Steps 3 and 5 ta^e 0(n^' tirre as they use pctr-s&l 

step 8 tal-es 0(n > tirre This can be done by the rrethod of 

parallel prefa> coirputation 

The forrration of sets X F I can be done in lines siirilar to 
the one used to fortr in the par-sel (as described in the 

analysis of step 4 of pars&l) and hence can be done in 0(n ) 

t lire 

We discard all together in steps 2 and 9 atleast one quarter 
of the points Therefore step 10 tales a tirre of t(3n/4> 
t(rO ~ t(dn/4) + c^n^ = 0(n^> 

We now address the question of recognising on what side of 

the straight line the unconstrained center lies The function 

2 a 

f(/ y) ~ rra {(> a^ > (y 1 < i S nl is convej- f is 

conve>' not only in each variable but also as a function of two 
var lab les 

h(y> = corve> 

By mnirrising gC ) we evaluate h(Ln The y coordinate of the 
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urcc ii&trained center is precisely where hCy) attains its 

mniiTuiT Since h(y) is convex we can tell the sign of y*" sirrply by 
looting in the neighbourhood of y = 0 Let ( >‘^ be the center 

in the constrained case Let 

I = { 1 C>* = g( 

If I = {il then > and y*" has the sign of If I = { i j I 

then (> 0) lies on the perpendicular bisector of the line 

c 

segrrent t (a b ) (a b )] y has the sign of the y coordinate 
11 0 

of the [Tidpoint of this segrrent i e l/2Cb^ + b^t In general all 

the points fa^ b^) i c I lie on a circle centered at ( ^ 01 If 

( (H IS in the conve> hull of these points then y*" = 0 

Otherwise there e> ist two points (a b , ^ (a b ) <i j c I) 

1 i J J 

such that f()- y) decreases as we rrove frorr ix in the 

direction of the mdpoint of the line segrrent C f a b '> (a b >1 

XX J 

(i e along the perpendicular bisector of that segrrent> In this 
c 

case y has the sign of (b + b )/2 The deterrr mat ion of these 

1 J 

E 1 

two points can be carried out in O^n ) tirre using O^n ) 

processors as follows 

Ai-ply affine tr ansf or rrat ions to transforrr ( and (a^ 

bj'> to <(i (H and hi 1) re=^pect ive ly If hi ii^ is not ir the 

convex hull of the given points then 3 an a such that ota^ 

for every i e I i e a rrust satisfy 

tTa> { b /a a iil a irin ih /a a, fi) 

ill 111 

and also b^ 0 for every i such that = (i If this is the 
case the points needed by us are the points at which the rra in utr 
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ar d the rriniirurr are obtained on the LHS and RHS respectively Note 
that this rrarirrutr and the rrinirrurr points can be obtained easily 
using 0(n^ processors in O(n^) tirre If (a^ b^) fa^ b^ ^ are 

the points so found then the sign of is that of l/2(b^ + 
Otherwise y^ = D 

Thus we can deterrrine on what side of the line the center of 

6 1 € 

the srrallest enclosing circle lies in 0(n ) tirre using OCn 1 
processors If the center lies on this line then we discover its 
exact location during the procedure 
Algor ithrr for the unconstrained case 

The following procedure salves the mnitrurr spanning circle problerr 
using 0(n^ CREW PRAM processors in Ofn^l tirre 
pcLr~uncons-pointsCSJ> 

(1^ If 1S| IS a constant nuirber use the sequential algor ithrr Else 

B 

subdivide the set of input points into subsequences of si e n 
and assign a processor to each subsequence 
For 1 = 1 to n^ ® do in parallel 

P^ pairs up the points in its associated subsequence and finds 

the bisectors 1 s of the line segrrents between the pairs of 

J 

points 
end for 

(2) Let L be the set of all the bisectors found above 

Divide the lines into two sets with the lines rra^ ing 

angles ot {3 with the > a is where ki < a < ||/2 and \\/2 < ft < 


2h 


0 



(3) Find the line of rredian slope (angle « with the > a is^ atrong 

rr 

the lines in the set using par-s&l Divide in to two 

sets one with slopes the rredian slope and the other with 

slopes the rredian slope Repeat the steps upto step 9 with 

the set L 2 tal- ing the rredian slope line of the set 
1 s 

(41 For 1 = 1 to n do in parallel 

pairs up the lines in its associated subsequence so that 

each pair has one line with nonpositive angle and the other 
with nonnegative angle considering the linear transfer rrat ion 
that tal-es the a> is to the line y = o» 

rr 

For each pair (1, 1 ) 

^ J 

If (1. 1 ) are parallel to the x a is y. is the rrean of the 

r J K J 

y coordinates Otherwise find point of their 

intersect ion 
end for 
end for 

(5> Find the rredian of all such y, s using par—sel 

' J 

(6^ Test on which side of the line y = y the center lies using 

rr 

the rrethod outlined in the earlier section 

(7) If the center lies on the line y = y we are done 

rr 

If the center lies underneath y = y^^ 

X = { } 

1 G 

For a = 1 to n do in parallel 

P considers the pairs <1 1 > in its associated set 

a 1 J 

If 1 II 1 and y > y and 1 lies above y = y then drop p, 

1 J 1 J rr 1 rr r 
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whtre 1^ IS the bisector of P^) and lies underneath 

1 

1 

If 1 1 are not parallel and y > y then add y to the 

1 j ij rr 1 j 

set X 
end for 

(8^ Find the tredian x of X 

rr 

(9) Deterrrine on what side of the line = > the center lies 

rr 

using the rrethod outlined in the previous section 

(ICn If the center lies to the left of this line (the renaming 

cases can be dealt simlarly> 

1 B 

For i = 1 to n do in para] lei 

P considers pairs <1, 1 ] such that y, > y and y, ^ y 

1 t'J ijrr lj rr 

If 1^ forrrs a nonpositive angle with the positive direction 
of the > a IS drop the point which lies southwest of it 
ffig 2311 
end for 

(111 Let S be the set of rerraining points pa.r-xincons-lin&CS J> 
233 Analysis 

Step 1 4 7 10 tal-e D(n^) tirre because each processor in 

these steps spends linear arrount of tirre on its input 

Step tales OCn^i tirre by following the analysis of step 4 
in pcxr~s&l 

Step 3 tales O(n^) tirre since the median can be found in 
0(n®l tine and the division of the set into two can be done as in 
step 4 of pax'—s&l in Ofn^ltiire 




I 

I 

::3 

— 


2 3 1 


Circled podnta are to be drsc 


Steps 5 and 8 ta^e tiire because they use par—sel 

Steps 6 and 9 tale 0(n^l txrre by the rrethod given earlier 
We drop atleast one si teenth of the points in each iteration 
Therefore step 11 tales t(15n/16) tirre 
t(ni = t(15n/16) + c^n® = 0(n®) 

2 4 Minirrurr spanning circle for a set of lines 
2 4 1 The ProbleiT 

The problen is to conpute the smallest radius circle which 

intersects all the given n lines in the plane The radius of such 

a circle is called the intersection radius A sequential linear 

tiire algorithrr for this problerr has been proposed in CBM9n] 

Let 1^ == a^N + b^y + c^ - H i - 1 2 n be the given set 

of lines in the plane The bisectors of their two lines 1 1 i!»! 

1 J 

j are given by 


(a + b y H c I 
1 1 1 


(a + b y + c ) 
J J J 


(a ■^+ b 


, 2^ , 2,0 5 
fa + b ) 

1 1 


We denote these bisectors by b and b 

1 J 1 J 

2 4 2 The Algorithrr 

We first consider the constrained case where the center of 

the circle is constrained to lie on a given line L We let this 

line to be the > ai is The following procedure conputes the center 

1 ““6 

of the rrinirrurr spanning circle for a set of lines using Otn ) 


processors The input to the procedure is a set of lines S and a 
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line L (here the a is^ The output is the center of the circle 
pcur-cons- 1 1 n&C 

(1) If |S| IS a constant use the sequential algorithrr else 

subdivide the input set of lines into subsequences of sice 

and assign a processor to each sequence 

1 6 

(21 For 1 = 1 to n do in parallel 

finds the intersection points of these lines in its 

cciated subsequence with the > a is 
end for 

(31 Let X be the set of the coordinates of the points found in 
step 2 above 

= par-selCX \X\y2:> 

(41 D - { } 

1 e 

For 1 = 1 to n do in parallel 

P finds the distance d =• distance of (> hi frorr L , for 

1 1 j fi J 

each L e the set of lines associated to P and adds d to D 
J 1 1 J 

end for 

(5) g = par-s&lCD \D\^ 

(61 I = i } 

1 6 

For 1 = 1 to n do in paral 1p1 

P adds J to I if dist(( Ui LI = g 
1 rr J 

end for 

(71 Initialise right - false left = false center = false 
If |I| n^ u e one processor 
else 
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Subdivide I into subsequences of size and assign each 

of these subsequences to a processor 

1 e 

For 1 = 1 to n do in parallel 

initialises right^ == false left^ - false centre^ = 
fa 1 se 


P sets right = 

1 1 

tr ue 

if every 1^ e 

its 

assoc lated 

subsequence touches 

the 

right sem-circle 

(Utv 

in 

fig 

241) centered at ( 

!i) 

and sets left = 

n 

tr ue 

if 

the y 


touch the left sem circle f usv in fig '^41) 
end for 

(81 Set left === true if all left true 

1 

Set right = true if all right^ = trie 
else return ( 1)) 

IT 

(9) F = { } Y = { } 

1 "~G 

For 1 = 1 to n do in parallel 

For each 1 in its associated subsequence P does 

if right = true then if 1 intersects a is to the left of 

J 

( (I) add 1 to F 

IT j 

if left =■ true then if 1 intersects a is to the right of 

J 

(> 0) add 1 to T 

rr J 

end for 
end for 

dtn Subdivide F into subsequences of si e and 

assign a processor to each sequence (Repeat the steps with 
the set T ) 
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For X 1 to ^ do in parallel 

P pairs up I s in F as , U > such that 1„ , 

intersects the > a is to the left of the intersection of 1^^ 

with the > a IS Let a be the intersection of 1 with the 

J J 

a is For each pair '^idd the bisector d who-^e 

intersection with the a> is doesn t lie between , and a^ 

•L— J J. J 

to the set D dig 2 4 2^ 
end for 


(Id Find d the rredian of the int er sect ions of d s with the 
rf J 

> a IS using par-s&l 


(121 Deterrrine which side of (d ff) the center (x 

IT 

required circle lies as in steps 7 and 8 


(tl of the 


(131 For 1 


1 to n 


1 e 


do in parallel 


For each d 


its associated subsequence does 


if X 
discard 1 
if 


d then if d intersects 
J 


a!' IS to the left of d. 


2j 1 


(fig 


rr 


rr 


then if d intersects a is to the right of d. 


discard 1^ (fig ^ 4 31 

^ vj 


end for 


(14) Let S = set of refraining lines pan — con£i-l-Ln»CS ;> 

Analysis 

let t(ni be the tirre for par-cons-line 
Step 1 tales constant tirre 

Steps 2 4 6 7 9 Id and 13 tale O(n^) tirre as each 
proce sor has to spend a linear arrount of tirre on its input for 
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Fig 2 4 1 


HG (X 






Pig 2 4 3 If X* < discard I 2 J 


I 


ate 


hpiae steps 


Steps 3 5 and 11 tale O(n^) tirre by the analysis of pcu — s&l 

Step 8 tales titre This can be done by the rrethod of 

arallel prefi cofrputatian 

step tales D(n itirre by the rrethod giver earlier 
Altogether n/S of the lines are discarded frotr consideration 
Therefore the tirre for this step is tf7n/8> 
t (r 1 = t (7n/8l + c^n® 

- O^n^l 

Unconstrained center c a^e 


I f ( D1 IS the center of the constrained problerr using 
all the processors in parallel find the set 


I € 1 

|a + c 1 

1 t 

2 

/ Ca ^ + 

1 

b 

1 

= g ( > ^ 0 ) } 




If I IS 

a singleton 

et then 

•o- 

y = 

a^ and y'" ha 

the 

Sign 

of 

b If I has 

1 

e actly two 

eleirents 

then 

y^ has the 

sign 

of 

the 


y coordinate of the intersection of the cor respond mg lines Else 
if the lines corresponding to the indices in I forrr a bounded 
cci ve _et enclo mg ( then this is the required center 


c 

cil erwisc y has the sign 

of 

the 

y-coordmate 

of 

the 

intersection point of the two 

1 mes 

such 

that to 

find 

the 

tr ue 

center we need to rrove towards 

the 

intersect ion 

point 

of 

the-=^e 

lines along the intersection 

point 

of 

these 

1 mes 

along 

t he 


bisector passing through (x'°' (.n 

The following procedure determnes the center in the case of 
the unconstrained center probletr 


r-iy 



>ar--%fncons!:-l in&CS^ 


1) If |R| IS a constant use the sequential algorithrr Else 
sjbdivide the input set of lines into subsequences of si e 
and assign a processor to each sequence 


C2^ For 1 = 1 to n 


1 e 


do in parallel 


finds the angles each of the lines in its a sociated set 
rial es with the > a is 
end for 

(3) Divide the lines into two sets L 


1 


with the lines in 


therr maling angles a ft respectively with the > a is where (f < 

a < [f/2 and ff/2 < ft < 

1 2 

(4> Find the lines 1 1 the lines of rredian slope arrong the 

rr rr 

Lj and u respectivdly 

(5) Rotate the > a> is so that 1 ^ is the new > -a is ^Repeat this 

rr 


and the following steps with 1^ too^ 

(6) Use all the processors to pair up the lines Cl, 


1 ^ so that 
J 


one has a ve and the other has a positive slope Find f 

y 1 the intersection point of these lines 
1 J 

(7> Find the rredian of s 

tr ij 


1 J 


> the center lies using 
rr 


(Si Te~>t on which side of the line 
the rrethod outlined above 

(9) As un c the center lies to the left of the line (The other 
ca'-es can be dealt simlarly > 

Find y the rredian of y s which lie to the right of the line 

IT ij 



iri) Determine the side of y = the true center lies A sues 
w 1 o g that it lies below this line ^The center lies in 
the quadrant LL (fig 2 4 4) ^ 
in B = { } 

!*• 1 "“G 

For 1 = 1 to n do in parallel 

coTisiders all pairs of lines (1 1 ) m its associated 

set whose intersection points lie in the quadrant UR and add 
therr to B 
end for 

Repeat the above proce<=-s with the set of lines in B which 
intersect LL to localise the center in one of the four 
quadrants created by the rredian lines Assurre that this 
quadrant is the lower left quadrant L L let U R be the 
upper right quadrant 

(ld» Use all the processors in parallel to discard one of the 
pairs of lines in the given pairs whose bisectors donot 
intersect 1 L pj L L (fig 2 4 5) 

(14) Let S be the set of rerraining lines paT-uncons-l n&CS J> 
243 Analysis 

^tep 1 tal-e<=" constant tirre 

Steps 2 5 6 11 12 and 13 tal.e □(n®) tiire 

step 3 tales 0(n ) tirre (by the analysis of step 4 of 
par-sol ) 

Steps B and Id tale O(n^) tirre by the rrethod given earlier 
The nurrber of lines rerraining for considerat ion are 
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3 




i27n/128 Therefore step 14 ta^es t(127n/n8^ 


t(n) - c^n + t(127ri/128i 
= OCn®) 


2 5 Minirrurr Spanning Circle for a set of line segrrents 
The problenr is to c.orrpute the sirallest radius 
intersects a given set of n line segrrents in the plane 
solved by reducing the problerr to that of corrputing the 
radius circle for a set of lines and points and use a 


This 


prune 


search technique [JM9(n The following subprobleirs are 


ived 


1 Given a set of n lines and points in the plane 
srrallest radius circle intersecting the e 


the 


2 Given a set of n rays in the plane show how these 


can 


converted to 

a set of n points and 

1 ines 




Then each line segrrent 

can be 

thought 

of 

as 

a ' 

intersect ion 

of two rays 

(fig 

2 5 1) 

and 

hence 

subprobleirs 1 

and 2 






We define 







g^ (> y > - 

n a { d 

( l<i<nj^ 

rra fd^ 


y)) 

1 

1 

in 

LI 

y> = 


( y ) ) 

Pi 

in 

PJ 


i l<i<n^^ 

4U 






g( y) = 

rra 


y )) 




where L is a 

set of lines P 

as a 

set of points 


(1 ^ 

the distance 

of the point ( 

y) frorr 1 and 

d^(p 

( 

y> 


cjggenen 


be 


ate 


we 


have 


■>) 


is 


is 

the 


distance between points p and ( y> 




Solving subproblLfT 1 


This probleiT is solved in sirrilar" lines as in the case of the 

problerr with points and lines dealt with in the earlier sections 

We pair up the points of P arbitrarily and corrpute their 

bisectors Let be the rredian of the intersection of the 

bisectors with the a is We corrpute g( P) the radius of the 

n 

srralle'-t circle centered at (y p) We determne the side of (> 

rr rr 

P) on which the true center lies depending on the points of 
contact of the lines which touch the circle and the points which 

lie on the circle Then we can throw away a fraction of the lines 

and points as outlined in the previous sections As given in the 
earlier sections we can deterrrine which side of the given line 
the true eerier lies The unconstrained problerr i_» salved in rruch 
the sarre way Each pruning step consists of throwing away about 

l/16th of the points and l/64th of the lines Since all this can 
be done in OCn^) tirre using 0 n^ processors by the results of 
the previous sections we can say that subproblerr 1 can be solved 
in O(n^) tirre using 0(n^ processors 
Solving subproblerr 2 

Consider the constrained case first For each ray r^ consider 

the line 1 noriral to it and passing through its tail ^fig 

1 

2^2) We conpute the rredian of the points of intersection of 
these normals with the s a is and locate on which side of this 
rredian the true center lies Suppose the true center lies to the 
right of this iredian Then for each of the nortrals which 
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nter'^ects the > a is to the left of this rredian point we replaice 

he corresponding ray by a line or a point If the ray and the 

edian point lie in the sarre half space of the two half spaces 

eterrrined by the norrral then the ray is replaced by its line of 

upport else by its tail point This can be done in 0(n®^ tiirs 
1 e 

ising O^n ) processors by allowing each of these processors to 
wort- on subsequences of rays feach of si e n®^ 

After doing this do the pruning step as indicated in the 
solution to the subprobleir 1 Iterate until a constant nurrber of 
lines and points retrain 

Tc solve the unconstrained case we consider the norrrals of 
the rays Rotate the > a/is to be parallel to the direction of the 
median slope The iredian slope is the line whose angle with the 
> a IS IS the rredian of the angles of the lines with the a>is 
This rotation is to be done twice The first tirre we consider all 
the lines which rrahe an angle a with the > aris where P S a < n/2 
The second tirre we consider all the lines which tral-e an angle (3 


with the N a 

IS where |f/2 

< ft < P Pair 

up 

the norirals 

With 

the 

ve and +ve 

slopes arbitrarily and 

consider 

the 

SB t 

of 

intersect ion_* 

‘ .. 


Let > y be 

IT IT 

the 

rr ed lans 

of 

y , s 

1 J 

and 

nj ® respect. vely 

Locate 

the center w 

r t 


and 

y 



Let the center lie in the lower left quadrant LL (fig 253) 

Of the set of norrral pairs which have their intersections in 
the upper right (UR) quadrant one of each does not intersect 
quadrant LL We have a replacerrent strategy for the corresponding 
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rays of these nortrals If the ray and LL lie in opposite half 
espacGS determned by the norrral then replace the nay by its tail 
point else replace it by its supporting line (fig ^ Thus m 

a single iteration we replace 1/8 th of the rays by lines or 

1 Bv 

points Note that all this can be done in Qfn ) tirre using n 
processors by allowing each processor to wor 1 on a subsequence o 


51 e n 

After replacing a fraction of the rays we do pruni 
(subprobleiT 1^ on the set of lines and points generated so 

^ ) t 1 IT S 

Iterating the e processes we can solve the probleir m 0(n 


using 0(n 


pr ocessor s 





Fig 2 5 4 


Ray r is replaced by the supporting line of 
r^ Ray replaced by its tail point 
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CHAPTER 3 


COMPUTING MINIMAL LENGTH LINE SEGMENT OF EXTERNAL VISIBILITY 
3 1 Introduction 

Thougt the notion of wea^ visibility has been well studied 
the notion of e ternal weal visibility hasn t been given truch 
attention It is shown in CBT891 that given a convey polygon P 

the (Tinirral length line segrrent frorr which P is weally e>ternally 
visible can be found in 0(nl tine The algorithn is based on the 
solution to a geonetric ninimsation probleir 

In this chapter we give a parallel algorithn to conpute the 
shortest line segnent fron which a given convey polygon i-s we ^ly 
e ternally visible in 0(n log ni tine using Ofn ) CHEW PRAM 
processors In section 3 2 certain definitions and lennas are 
given In section 3 3 the parallel algorithn is proposed which is 
analysed in section 3 4 

3 2 Preliminaries 

For any integer n2:3 we define a polygon in the Euclidean 

plane a the figure P = Cp^ p 2 formed by n points p^ 

in E and n line segnents Cp^ ‘^i + l^ i = 1 2 n 1 

and Cp^ Pj^I The points p^^ are called the v&rtic^s of the polygon 

and the line segments are termed its &dges 

We assume that the vertices of P are in ^©nez'al position i e 
no three vertices are collinear and that the polygon is in 
standard form i e the vertices appear in counter cloct wise order 



A sirrple polygon P is said to be uieoMly externally visible 
frorr a line segirent L if L is outside P and for every point > on 
the boundary of P there is a point y on L such that the interior 
of the line segrrent Cr y] doesn t intersect P 

HL(> y) denotes the closed half-plane to the left of the 
directed line determned by the two ordered points The 
corresponding closed half plane to the right of the directed line 
deterrrined by two ordered points is denoted by HR(> yl 

The edge tp p ,3 of P is denoted by e If the half plane 
cct tains the interior of P it will be referred to as an interior 
half plane We denote the infinite half ray starting at point r 
and traversing a second point y by rayCx y3> 

The interior cone of support at p^ denoted by i-coneCp^^ is 
the wedge determned by HL(p^ ^ n HLCp^ + <fig 3 n 

A line of support of a corve> polygon P is a line that has at 
lea'*t one point in cotrrron with P but with the property that all of 
P lies on one side of the line 

A pair of vertices of the polygon that adirits parallel 

supporting lines will be called antipodal 

In fig 32 P^ ^nd p^ are not antipodal since parallel lines 

of support can t pass through p^ and p^ sirr ultaneously p^ and p^ 

fortr an antipodal pair 
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) Lone, 




Pig 3 



Illustrating the internal cone of support of 
vertex p of a convex polygon P 





Letrira 3 2 1 l'^'CP) is tangent to P 

Lerrma 3 2 2 L’^CPl has one of its endpoints on one of the 

bounding rays of i cone(p^) and the other on the other bounding 
ray for sorre value of i ^1 ^ i i n^ 

Lemfra 3 2 3 i P is weaf ly e ternally visible frorr LCa bl if and 
orl> if 3 tangent rays of support to P frorr a and b rayCa P> and 
ray(t P) such that the following three conditions hold 

(i) rayCa P) and ray(b P) intersect at soire point 

(ii) >■ IS a vertes of P and 


<iiii P IS contained in Aabx 

Lerrrra 3 2 4 A line segrrent L = [a bl lying in i cone(p^) 
tangent to P from which P is weaUy e ternally visible rrust tae 


tangent to a verte p^ such that p^ is an antipodal verte of 
Proof If L does not intersect ray (p^ P^ then edge Cp^ ^ 

P ] of P IS not visible from L Simlarly if L does not intersect 
the raytP^ P^^l> then edge Cp^ P. + i^ of P is not visible fro^ L 
Therefore U mu'-t have one end point on the ray Cp^^ P^ ' snd the 


other on the 
vertex of p^ 
tangent line 


rayCp Pj+i' assutre that p^ is not an antipodal 

Let P . s denote the vertex of P determned by a 
L-ei P^ (1-1 1 ) 

parallel to the rayCp^^ P^^ and antipodal to p^_j 


andP, S.n.Urly let P, ,, the vert, of P deterp.n.d 

by e tangent line parallel to raylPj P,^i' and antipodal to P^and 



Pi 4 l Clearly as we rotate a line of support in a 
cour ter cloc^ wise rranner starting at and ending at 

^Ki i4n polygonal chain C ^ ^ ^ ^ 

f 1 1 + 1^ 1 1 1 + 1)^ which has the property that all its 

vertices are antipodal to p^ Fu Ihernore if p is not antipodal 
tc p^ ther L rrust have an unoriented direction that lies in the 
wedge deteririned by the internal angle of P at the verte p This 
irrplies that one endpoint of L irust lie in int(i cone(p^)l whiuh 
IS a contradiction 

Leirrra 3 2 5 Let p be any point in the interior of the i cone of 
a verte p^ the length of the line segrrent CH K1 such that H and 
K lie on the boundary rays of i coneCp^l and passing through point 
P IS a unirrodal function and the rrinirrurr can be found in constant 
tirrc CBTSPI 

Lenrra 3 2 6 The total nurrber of antipodal pairs is N The set 
of all antipodal pairs can be found in linear tirre CSM753 
Proof The set of all antipodal pairs can be determned as 
follow Tran late the edges of the polygon to the origin treating 
theiT as vectors In this mapping edges go to vectors and vertices 
go to sectors (fig 3 31 In order to find the antipodal pair 
corresponding to sorre direction deterrrined by a line L translate 
the line so that it passes through the origin of the vector 
diagrarr The sectors through which it passes through indicate the 
points of the antipodal pair Deteririning the sectors which the 
line passes through can be done in 0(log n) tirre The set of all 
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antipodal pairs can be found out in 0(n) tiire by scanning 
sequettially around the vector diagrarr 


3 3 The Algorithrr 

Leirrras 321 322 323 and 3 2 4 irrply that any (a 

line frorr which a polygon is weaHy e ternally visible! rrust be 

tangent to a verte> p^ and have its end points on the bounding 

ray cf i core(p ! such that p and p forrr an antipodal pair 

+ 

Ir crdtr to find L (P ) through a verte p^ we need to 1 now 

the bounding rays on which its endpoints lie in order to apply 

lenna These bounding rays are deternned by i cone(p ) where p 

J J 

and p^ are antipodal pairs Our algorithm finds out all si h 

H 

po itle L tP! in parallel and finds the one of ^rinirrurr length 
arrong tho'"e 

The following procedures run on a shared aetrory SIMD rrachine 
with ^ proces ors P. 1 e where del 

poT-short-vts-ltno ctrifutcs the shortest line segrrent frorr which a 

1 e 

9 1 vf n conve> polyson is weakly e>ternally visible using n 
CREW PRAM proce <=^or 5 The input to the procedure is the polygon P 
(a ‘~et of vertices V and a set of edges E! The output is the 
shorte t lint segnent frorr which the polygon is weal-ly e ternally 
visible 

pcur—chort-vts-lmoCV 

(1! Subdivide the sequence of edges E into n subsequences each 
of Bi e n^ and assign a subsequence to each processor 



1 o 

( > For 1 1 to n do in parallel 

tran lates each of the edges in its associated subsequence 
tc the origin treating therr as vectors 
end for 

1 ““G 

(31 For 1 = 1 to n do in parallel 

deterfTines the antipodal pairs corresponding to the 
direction of each of the edges in its associated subsequence 
by determning the sectors through which the lines in those 
directions pa<^s through 
Cl d for 

(41 Fir k := 1 to n^ ^ do in parallel 

For eacti antipodal pair (p^ P^l associated to 

find the shortest line segrrent through p^ with its end 

points on the bounding rays of i coneCp i If this line 

J 

IS not a line of support to P (i e it cuts the interior 

of pi P^ tales the shorter of the line segrrents through 

(p , P 1 and (p P ,.l with the end points on the 
bcunding lines of x coneCp^i (fig 3 4i Pj^ simlarly 
find the line segnent through p^ P|^ then finds the 
shortest of the line segments that it has found using the 
antipodal pairs allocated to it and writes it in SCt-l 
the ^ PO->ition of an array S in shared rreirory 
end for 
end for 

(5) Call par-SfhorUS:> to find the shortest line segrrent 




The following procedure finds the shortest of the given set of 
line The input is a set of line segrrents along with their 
lengths 

The output 1 -. the shortest line segrrent 
pcxi — sf\or t C S-5 

1 If 1*^1 3 then u c one processor to return the shortest line 

else 

1 p 

Bui divide B in to |B| sub equences each of which 
ha 1S(^ elerrents and assign a subsequence to each 
pr ocessor 

For 1 = 1 to |S|^ ® do in parallel 

P>ttes or find-> the shortest s^ of its a sociated 

ut equerce u ing the sequential algorithrr for finding 1 
shor test 

P writes s in M(il the i^*^ position of an array M in the 
1 1 

shared rrerrary 
3 poo — shor tCMy 

3 4 Analysis 
Analysis of par-short 

Let tfnl be the running ti^e of par-short for an input of 
SI e n 

Step 1 tales constant tirre 

The sequential algorithrr to select the shortest is of linear 
tirre Therefore step 2 tales tirre units where Cj is a 


const ant 



Since |M| = n step 3 requires t(n^~®) time 

From the above we have 

t(n) = c^n^ + t(n^ whose solution is t(n) - 0(n®) 

Analysis of par short-vis-lrne 

Step 1 needs constant time 

6 

step 2 tales c^n time units where is a constant because 
tie translation of a single edge tales constant time and there are 
n® edges associated with each processor 

In step 3 the antipodal pairs are found by initially 
determining the antipodal pair corresponding to the direction of 
the beginning edge in the associated set of edges of a processor 
The antipodal pairs corresponding to the directions of the rest of 
the edges is found by a linear scan To find the antipodal pair 
corresponding to the first direction it takes Odog n) time 
according to lemma 3 2 6 To find the rest of the antipodal pairs 
it tal<es O(n^) time Hence the time for this step is dominated by 
n® Therefore the time for this step is c^n® where is a 

constant 

In step 4 each processor tahes constant time to find the 
line segment corresponding to one antipodal pair according to 
lemma 3 2 5 Therefore this step requires CgO® time units where 

c„ is a constant 

w 

, e ftma iintr*; where c- is a constant 

Step 5 needs c^n time unix-s wnere 

according to the analysis done above 
Therefore t(n) = 0(n®) 



CHAPTER 4 


RECONSTRUCTION OF AN ORTHOGONAL POLYGON FROM ITS 
VISIBILITY GRAPH 


4 1 Introduction 

The recogrition problerr for visibility graphs is given a 
grat h to determine whether it is the visibility graph of a sirrple 
polygon The conple ity of this problerr is un^nQwn In COJ873 a 
linear tirrp algorithrr has been given for the r econstr uc 1 1 on of an 
orthogcnal polygon giver the circular enbedding of its visibility 
graph In CF( 9U] a linear tirre algorithrr has been given for the 
recogi ition of the visibility graph of a spiral polygon 

In the present chapter a parallel algorithir is proposed for 
the problerr of constructing an orthogonal polygon given its 
visibility graph The algorithrr runs on Dtp) CREW PRAM pro essors 
in O((log p)ri/pl tirre In section 4 2 prel ir mar les are given A 
sequential algorithrr is briefly reviewed in section 4 3 In 
section 4 A a parallel algorithir is proposed which is analysed in 

section 4 5 

4 2 Prel iminar less 

In this section certain definitions and notations are given 
A polygon is a finite figure in the plane that is bounded by 

a finite nurrber of straight line segrrents 

A singly connected polygon is bounded by n points 



V called vc^rtic^s) and n line segrrents Cv, Cv^ 

j ■^rid Cv^ 3 (called &dges) 

Ai orthogonal polygon is a polygon whose edges are all 

aliyitd wilfi a pair of orthogonal co ordinate a es which we tale 

to tie the horizontal and the vertical without loss of generality 

Thus tt e edges alternate between hori ontal and vertical and 

always (rcet orthogonally with internal angles of either 90 or 

^/(i (Fig 4 1 IS an orthogoral polygon ) 

An orthogonal polygon is said to be in ^©neraZ. position if 

no two verti es can be connected by an internal hor i ontal or 

vertical lire segnent that doesn t intersect P s boundary 

Throughout the discussion here we use polygon to trean 

orthogcnal polygon in general position 

Two vertices v and v of a polygon P are visible if the 
X J 

clct td line segner t between therr doesn t intersect the e tenor of 

P Note that if the line segrrent touches the boundary of P the 

vertices are still considered visible 

Two edges e and e are said to be horiBontally Cvertically:> 
X J 

visible iff they can see one another along a vertical (hori ontaH 
line that is iff there exists a vertical (hori ontal^ line 
segrrent interior to P with end points on e^ and and which 

doesn t otherwise intersect the boundary of P dn fig 4 1 edges 
a and 4 are visible while edges 0 and T are invisible ) 

An edge is said to be a edge if the interior of 

the polygon lies aboveCbelow^ it Sirrilarly it is said to be a 
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lv/tCrtght!> edjc if the interiot' of the polygon lies to its 

l»/tCrtght^ (In fig 4 1 edges H and 4 are bottorr and top edges 

r ppctivply 

A vertical edge C(>j y^) (k^ y^) 1 is said to be go rxg up 

if y /o else it is said to be going down A hor i ontal edge 
C( j y^l <>2 said to be going left if i ^ ^2 

said to be going right 

The nodes of a visibility grap^ correspond to georretric 

cofTPonents ‘^uth a vertires or edges and two nodes are conne ted 

by an arc of tht graph if the cafrponents can see one another 

perhaps under sone restricted fortr of visibility 

The orthogonal edge visibility graph. G for an orthogonal 

polygon P is defined as follows 6 contains a node for each edge 

of P and twe nodes a&->Dciated with hori ontal (vertical l edges e^ 

and e are connected by an undirected arc in G iff they are 
J 

hcii crtallyCverticallyi visible In the discussion here the tern 
visibility srsih IS used for the terr orthcorel edge 


visibility graph 

Tt-e visibility graph of an orthogonal polygon con ists of two 
disjoint trees one cor respond ing to the hori ontal visibility 
(called the horiisontal visibility tree and the oth 


corresponding to the vertical visibility tree y 

we say that a tree is realisable if there is a polygon with 
the tree as one of the two corrponents of its visibility gra^ h 

We say that two trees can if they are jointly realisable 



by the 


ti polyg(3n 


A lat>t^ltne cf a tree of n nodes is a bijection between th 
nc tif ar d the set of integers 0 1 n-1 

A labeling of a tree is realisable if there is a polygon that 
•'t. ill f tie tree and such that the polygon edges [ray be nurrbered 
f'* 1 ^1 ^ counter cloclwise traversal of the boundary to 

agree witli the labeling 

Defire d circl© embedding of a tree as a layout of the tree 
witl ir a circ le in the plane such that each arc is trappped to a 
chord of the circle and all nodes of the tree are trapped to lie 
on tlic circle (fig ^ '^1 The circle corresponds to the boundary of 
the pilyjrn inflated to a circle and an etrbedding is 
tor r Ic jic al 1 y equivalent to a layout of the visibility trees 
within the polygon it elf A labeling of a tree maps to a unique 
errbedd i ng 

The portion of the polygon that reali es the subtree from a 
node to its iirirediate descendants is a histogr-am or Manhattan 

slyline polygon This polygon is the orthogonal edge visibility i 

r 

polygon for edge e enclosing all the points visible to edge e 

Two nodes of an errbedded tree are said to be 2-adj<xc»rti or r 
ixf&xghhoxirs if they are adjacent on the circle if the other tree is I 
ignored Ncte that 2 adjacent nodes receive labels i and i+2(rrod 

n) in the labeling I 

The distance between two 2 adjacent nodes is defined as the 
nufTber of arcs in the path in the tree between them Let be the 

45 I 




di<=itaTicG betwcer node i and its ccunter cloc^wl 5 e neighbour 


Lerrrra 4 2 1 In a reali able embedding the distance d^ between 
each pair of 2 adjacent nodes satisfies 5 3 Moreover the two 
angles in a realizing polygon between two adjacent nodes are 
deterrrired by d^ follows where c and r rrear conve and refle 
angles respectively 

d = 1 cc 

1 

d =2 rc or cr 

1 

d = 3 rr 

1 

Px'oof - Let an arbitrary hori orital edge of a polygon which we 
talf to be a botton edge without los of generality be labeled n 
and label the reiraining edges with increasing mde 

counter -clod wi e The right end point of 0 is either a convex or 
a refle?’ verte? If this end point is conve then distinguish two 
further cases depending on whether the upper end point of 1 is 
convex or refle> The forrrer case is illustrated in fig 4 3<a> 
and justifies the clairred cor respondence between d^ = 1 and cc In 
the later case (fig 4 3(bd there rrust be an edge a above 1 

leadin 9 to ^ ^ind angles cr 

If the right end point of d is refle we again have two 

cases depending on whether the lower end point of 1 
reflex In the forrrer case (fis 4 3 (cn there is again an edge a 
above 1 and = 2 with angles rc In the latter case (fig 

4 3(d) ^ there rrust be an edge a above 1 and an edge b 
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c 



( a) 


(b) 



h 

(c) 


c 



d 

(d) 


Fig 4 4 


Illustrating projection constraints 



This ■algorithn' identifies which edges of the orthogonal 
polygon determne the vertical edge- of the edge visibility 

polygon Once these are identified we 1 now their left to right 

ordering 

Let e be a bottorr edge whose node ha degree greater than 
one and let a and b be two edges visible froir e with the edges 
occuring in counterclocl wise order e a b and with no edges 

between a and b also visible to e Then the vertical edge of e s 
edge visibility polygon between a and b is detertrined by either 
a+Hrrod n) or b Hrrod n) Which case obtains can be decided by 

checi’ ing whether node a+Hnod n) is connected by an arc to any 
node in the range b countercloct wise to e If node a+Krrod n) is 
connected to a node in the range b counterc loci- wise to e then 
a+ltrrod nl is the vertical edge of e s edge visibility polygon 
else b-l(rT)od n) is the vertical edge Continuing in this rranner 
we can conclude that e,\ e„<. e where a ■=' b treans that the 

1 ^ j 

> coordinate of edge a is less than that of b and e. e^ e 

I ^ J 

are the vertical edges of the edge visibility polygon of edge e 
where e is a bottorr edge If e is a top edge the ordering would 
have been e^^ e^ e^ After all these partial orders have been 

delernined they are rrerged The position of an edge in this list 
IS used as the coordinate for the vertical edge Sirrilarly the 
y coordinate for the hori ontal edge can be deterrrined Thus the 
polygon can be ronstrut 'v-’d 


I 
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4 4 The parallel algorithrr 

A -.traight forward paralleli ation of the above (tethod would 
not yield an efficient algorithfT Instead we ta^ e up a different 
approach We detertrine the angle sequences of the polygon and the 
length of the edges frofr the given circular eirbedding of the two 
trees Letrras 4 2 1 and A 2 2 are used to determne the angle 
sequences and the edge lengths are ta^en to be the degrees of the 
corresponding nodes in the graph 

The following algorithrr runs on a shared uerrory SIMD rrachine 
witti F processors F. P„ P where 1 p n The input to the 

1 P 

aljorithrr is the circle enbedding of the visibility graph tan 
array of adjacency lists^ The coordinates of the vertices of the 
orthogonal polygon (whose visibility graph is the input) 
constitute the output 
par~vis~real 10 

(1^ Subdivide the input sequence into p subsequences each of si e 
10 length of the input list) / p and assign a prouess»ar to ea h 
subsequence 

(2) For 1 = 1 to p do in pavallel 

P scans it associated subsequence for the beginning of the 
1 

adjacency list of ary elerrent and if it finds one it stores 
the corresponding array inde in the shared rreirory so that 
further acccs_t to the adjacency lists of the elerrents can bd 
done in constant tirre It also stores the elerrent to whos^ 
adjacency list the elerrent in the beginning of its associated 
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subsequence belongs to 
end for 

(3) For 1 » 1 to p do in parallel 

For^ ton/p do 

finds the degree of the node (i-1)n/p+|' by using the array 
indices found in step 2 stores the degree of node 

(i-1)n/p+h as the length of the edge (i-1>n/p+h in the shared 
IT err 0 ry 
end for 
end for 

(4) For 1 = 1 to p do in parallel 

For h = 1 to n/p do 

P^ initialises dC ( l-1 ) n/ p+h H = 3 and 

pathC( 1 - 1 )n/p4h 3 = 0 

end for 
end for 

(5) For 1 = 1 to p do in parallel 

P^ processes the subsequence associated with it 

While scanning the adjacency list of node a if it encounters 
a node b such that b - a + 2 (trod n) then it sets dCa3 == 1 
If it encounter nodes b and b + E(n'od n) in the adjacency 
list of node a then it sets dCb3 = 2 and pathCbl = a 
end for 

(b) For 1 = 1 to p do in parallel 

(6 1) P^ determines the angle sequence of the polygonal chain 
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frorr edge (i“ 1 )n/p +1 to iCn/p) as follows 
For each l< in this range 

If dllh3 = 1 then the angle sequence is cc 

If dCh3 = 3 then the angle sequence is rr 

else if node h + 1 projects across edge path(k)) 

then the angle sequence is cr 
else the angle sequence is rc 

(6 2) determines whether edges (i-1)n/p+S to i(n/p) are 
top bottom left or right assuming edge (i-l)n/p + 1 to 
be a bottom or a left edge and using the angle sequence 
found above 

(6 3) determines the coordinates of the vertices of the 
edges (i- 1 )n/p +1 to i(n/p) ta^ mg the coordinates of the 
beginning verte< in the chain of edges to be <0 0 ) 

using the angle sequences found in step 6 1 and the 

lengths of the edges found in step 3 
end for 

(7) Assign n/p edges to each processor 

For 1 = 0 to log p - 1 do 

■ *1 ^ *| 

For [ = 0 to n - 2 in steps of 2 do in parallel 

P tests the beginning edge e, in its associated set of 

h+ 2 ^+ 1 

edges and the la t edge e ^ in the set of edges associated 

to processor P Let e. be a vertical edge and e^^ be a, 

h+2^ 

hori ontal bottom edge going right 

Let = C( ^ y^) ( ^ = C( ^ yg) (<3 y 3 ) 3 
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Case* 


I 16 a left edge going up and the angle between the 

edges is conve A nrirror irrage w r t y = is to be 

ta^en (fig 4 5a) 

II e^ IS a left edge going up and the angle between e^ and 

I'efle A mrror irrage w r t ) ^ is to be 

ta( en (fig 4 5b> 

III e^ IS a left edge going down and the angle between the 
edges is refle A rrirror irrage w r t y = y^ and with 
respect to = > ^ is to be tahen (fig 4 5c) 

IV Cj IS a left edge going down and the angle between the 

edges is convei No rrirror irrage need be tah en 

Similarly the othrer cases can be dealt with 
1 1 i 

For j «= h + 2 + 1 to I + 2 do in parallel 

P translates the coordinates of the vertices of the edges 
in its associated sequence Then it tahes the rrirror images 
of the vertices if needed as determned above 
end for 
end for 
end for 

4 5 The Analysis* 

Step 1 tah es constant tirre 

Steps 2 and 5 tah e c^n-'p and c^n/p tirre respectively where and 


c„ are constants 

iL 





step 3 lal- es C 2 (n/p) litre 
Steps 4 and 6 la^ e c^n/p and c^n/p tiir 
and are constants (In step 6 1 to 
across edge (i patl-(i)) it ta^ es cor 
4 2 £ ) 

Step 7 ta^ es c,(log p)n/p litre where c 
o 

processor tahes a constant litre to tr 
nrirror iirage of the coordinates of one v 
t(n) = c^n/p + CgO/p + c^n/p + ^ 4 ri'^P 
= 0( (log p)n/p) 


S3 



In the above algorithm taking p = n / log n we get an optimal 
parallel algorithm as follows In step 7 we can do the merging by 
doing a prefix computation of the amount of translation to be 

done ( 5 imi la r ly the necessity of taking mirror images) Hence this 
step can be done in 0(lag n) time using 0(n / log n) processors 

Therefore the time complexity of the algorithrr is O(log n) 

Another method to determine the coordinates of the vertices 

of the polygon is as follows The direction in which an edge goes 

changes when two convex or reflex angles occur together We note 

the occurence of such a sequence by -1 and the occurences of 

other angle sequences by +1 and then do a prefix multiplication on 

these Then multiplying the numbers so obtained by the edge 


lengths 

we get 

the 

relative 

lengths of 

the edges Doing a 

pref IX 

sum computation 

on 

these 

gives the 

exact coodinates 

of the 

vert 1 ces 

Using 

0(n 

/ log n) 

processors 

this can be done in 

Odog 


n) time 



I 

CHAPTER 5 

CONCLUSIONS 

In this thesis parallel algorithms for certain problerrs ] 
Con |i ut at lonal Gecrretry have been presented The shortest lit 
segrrent from which a given conve> polygon is weafly external 
visible is coiTputed in 0(n^ tirre using 0(n^ CREW PR/ 

processors The reali ation of the orthogonal polygon given tt 
circular embedding of the two visibility trees of the polygon c 
be done in 0((.n/pl log titre using CREW PRAM processor 

The linear prograrrir ing problerr the rrmirruir spanning circ 
problem of a set of points lines and a set of line segments ha 
been solved in 0(n®> time using Ofn^ CREW PRAM processors 

I 

Several problems remain open One is the problem of findi 
an efficient parallel algorithm for the computation of tl 
shortest line segrrent from which a *^imple polygon is wea^ 
externally visible Another is to find an optimal parall 
algorithm for the r econstr uct ion problt of a simple polygon giv 
its visibility graph 
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